Counterflow of spontaneous mass currents in trapped spin-orbit coupled Fermi gases 
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We use the Bogoliubov-de Gennes formalism and study the ground-state phases of trapped spin-orbit coupled 
Fermi gases in two dimensions. Our main finding is that the presence of a symmetric (Rashba-type) spin-orbit 
coupling spontaneously induces counterflowing mass currents in the vicinity of the trap edge, i.e. f and 4- 
particles circulate in opposite directions with equal speed. These currents flow even in noninteracting systems, 
but their strength decreases toward the molecular BEC limit, which can be achieved either by increasing the spin- 
orbit coupling or the interaction strength. These currents are also quite robust against the effects of asymmetric 
spin-orbit couplings in x and y directions, gradually reducing to zero as the spin-orbit coupling becomes one 
dimensional. We compare our results with those of chiral p-wave superfluids/superconductors. 

PACS numbers: 05.30.Fk, 03.75.Ss, 03.75.Hh 



I. INTRODUCTION 

The controversy over the expectation value of the intrin- 
sic ground-state angular momentum of the A-phase of su- 
perfluid 3 He in a given container has a long history, and 
it has not yet attained a universally agreed resolution |Q]]. 
For N weakly-interacting particles, theoretical predictions 
for the BCS ground state vary orders of magnitude from 
N\A\h/(2e F ) to Nh/2, where |A| is the amplitude of the 
order parameter, i.e. pairing energy, and ep is the Fermi 
energy. At a first glance, the former expectation, which is 
based on the observation that pairing affects not all but only 
a small fraction (\A\/ep)N of fermions, seems more in- 
tuitive. However, this subject has recently regained inter- 
est in condensed-matter and atomic and molecular physics 
communities due to its relevance to the possible p x + ip y - 
superconducting phase of Si"2Ru04 [2] and atomic chiral p- 
wave superfluids [3], respectively, supporting for the latter ex- 
pectation (j, |5J . On one hand, these recent results seem intu- 
itive at least in the strong fermion attraction limit of tightly- 
bound p-wave Cooper molecules, where all of the molecules 
can undergo BEC at zero temperature with each molecule hav- 
ing a microscopic angular momentum h. On the other hand, 
they seem counterintuitive in the weak fermion attraction 
limit of loosely bound and largely overlapping p-wave Cooper 
pairs, since these results imply that the angular momentum 
in the BCS ground state is macroscopically different from its 
vanishing value in the normal ground state, when the pair- 
ing energy becomes arbitrarily small. While this controversy 
still awaits for an experimental resolution, here we propose an 
alternative atomic system where the microscopic angular mo- 
mentum of Cooper pairs can also give rise to a macroscopic 
one. With their unique advantages over the condensed-matter 
systems, it is plausible that the macroscopic angular momen- 
tum of this alternative system may be observed for the first 
time with cold atoms, and that this would also shed some 
light on the 3 He controversy for which the basic mechanism 
is found to be similar. 

Cold atom systems have already proved to be versatile in 
simulating various many-body problems. For instance, one 



of the major achievements with atomic Fermi gases in the 
past decade is the unprecedented realization of the BCS-BEC 
crossover where the ground state of the system can be 
continuously tuned from the BCS to the BEC limit as a func- 
tion of increasing fermion attraction with the turn of a knob. 
The main difference between the BCS-BEC crossover and the 
BCS theory is that the Cooper pairing is allowed not only for 
fermions with energies close to the Fermi energy but also for 
all pairs with appropriate momenta lil|6|]. Having established 
the basic tools, it is arguable that one of the most promising 
research directions to pursue with cold Fermi gases is the ar- 
tificial spin-orbit gauge fields. Such a field has already been 
realized with cold Bose gases using a novel technique 

based on spatially-dependent optical coupling between inter- 
nal states of the atoms, and the same technique is equally ap- 
plicable for neutral fermionic atoms. In fact, we have recently 
learned that the first generation of quantum degenerate spin- 
orbit coupled Fermi gases has recently been created with 40 K 
atoms H]. 

This immediate possibility of creating a spin-orbit cou- 
pled Fermi gas has already received tremendous theoretical 
interest in condensed-matter and atomic physics communi- 
ties, exploring a number of exotic superfluid phases with bal- 
anced or imbalanced populations, at zero or finite tempera- 
tures, in two or three dimensions, etc. lllOl - IZJ . Motivated by 
these recent advances, here we study ground-state and finite- 
temperature phases of trapped spin-orbit coupled Fermi gases, 
and our main results are as follows. We find that the presence 
of a spin-orbit coupling (SOC) spontaneously induces coun- 
terflowing t and X mass currents in the vicinity of the trap 
edge 12411 . We show that these currents flow even in nonin- 
teracting systems, and that they are quite robust against the 
effects of asymmetric SOC in x and y directions, gradually 
reducing to zero as the SOC becomes one dimensional. How- 
ever, their strength decreases toward the molecular BEC limit, 
which can be achieved either by increasing the SOC or the in- 
teraction strength. We argue that the origin of spontaneously 
induced counterflowing mass currents in spin-orbit coupled 
Fermi gases can be understood via a direct correspondence 
with the chiral p-wave superfluids/superconductors, however, 
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with some fundamental differences in between. 

The rest of the paper is organized as follows. In Sec. [II] 
first we introduce the mean-field Hamiltonian that is used to 
describe the spin-orbit coupled Fermi gases on optical lattices, 
and then use the Bogoliubov-de Gennes (BdG) formalism to 
obtain the generalized self-consistency equations for the num- 
ber of fermions, superfluid order parameter, and strength of 
circulating mass currents. The numerical solutions of this for- 
malism are presented in Sec. [TTTJ which is followed by a brief 
discussion on their experimental realization in Sec. [IV] We 
summarize our main findings inlVl 



II. THEORETICAL BACKGROUND 

In this paper, we consider a harmonically trapped two- 
dimensional optical lattice in the presence of a non-Abelian 
gauge potential. In addition to the usual spin-conserving hop- 
ping terms, the main effect of such a gauge potential is the 
appearance of an additional spin-flipping hopping term in the 
single-particle kinetic energy term as discussed next. 



A. Hamiltonian 

Spin-orbit coupled Fermi gases on optical lattices can be 
described by the grand-canonical mean-field Hamiltonian, 

H = ~ *E ( C l+e^+e^ + H.cZj 

i,e 

- J2 (/v + v£ - vf ) c\ lCai 

i 

+ E ( A 44* + A *wi) . (d 

i 

where the operator <y ai (c„i) creates (annihilates) a pseudo- 
spin a = {t, 1} fermion at lattice site i, the spinor C\ — 
( c |i' c L) denotes the fermion operators collectively, e = 
{x, y} allows only nearest-neighbor hopping, and H.c. is 
the Hermitian conjugate. For a generic gauge field A = 
(aa yi —f3a x ), where a e is the Pauli matrix and {a, /?} > are 
independent parameters characterizing both the strength and 
the symmetry of the spin-orbit coupling, the f and J. fermions 
gain <pi+x,i = e^ lo " 7y phase factors for hopping in the posi- 
tive x direction and 4>i+y/ L — e i,9crx phase factors for hopping 
in the positive y direction. Note that the spin-conserving and 
spin-flipping hopping terms are proportional, respectively, to 
cos a and sin a in the x direction, and to cos /3 and sin /3 in 
the y direction. The off-diagonal coupling Aj is the mean- 
field superfluid order parameter to be specified below, and 
Vjf = gri- a i is the Hartree term where g > is the strength 
of the onsite attractive interaction and n a i is the filling of a 
fermions at site i. Here, (— y) —i an d vice versa. We in- 
troduce er-dependent chemical potentials /u CT to fix the num- 
ber of a fermions independently, but we assume both t and 
I fermions feel the same trapping potential V t T = Vbr?/2, 



where the distance r; is measured from the center in our Lx L 
lattice. Next, we solve this Hamiltonian via BdG formalism. 



B. Bogoliubov-de Gennes formalism 

The BdG equations are obtained by diagonalizing the 
quadratic Hamiltonian given above via a generalized 
Bogoliubov-Valatin transformation, leading to a 4L 2 x 4X 2 
matrix-eigenvalue problem, 
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where u a ni and are the components of the nth quasi- 
particle wave function at site i, and e n > is the corre- 
sponding energy eigenvalue. While the off-diagonal cou- 
plings are Ay = A^y diagonal in site indices since we 
consider onsite interactions only, the nearest-neighbor hop- 
ping and onsite energy terms can be written compactly as 
rj a , = - ( Mff + V* - ^ T )%<W, where % is the 

Kronecker delta. Here, the non-vanishing nearest-neighbor 
hopping elements are t l ^ x — t cos a and t l ^ +x = — t l ^ +x = 
—t sin a for the positive x direction, and t % ^^ — t cos (3 and 
i^'! +y = t l ^ +v = it sin P for the positive y direction. Note 
that the hopping in the negative x and y directions are simply 
the Hermitian conjugates. 

Equation (0 needs to be solved simultaneously with the or- 
der parameter A, = g{c^iC i% ) and number n al = (c^c^) 
equations for a self-consistent set of Aj and fi a solutions. 
Here, (• • • } is a thermal average. Using the Bogoliubov- 
Valatin transformation, we obtain 



Ai = <?E l^iiTfi-^ + uUvlXf^n) 
n 

n«i = E [K! 2 /(£n) + K^fi-Sn) 



(3) 

(4) 



where f(x) — \/{e x / T + 1) is the Fermi-Dirac distribution 
function and T is the temperature. Here, we set the Boltzmann 
constant ks to unity. Note that < n a ; L < 1 is the number of 
a fermions at site i (number filling) in such a way that N a — 
J2i n ai gives the total number of a fermions. Equations ([2]-|4]) 
correspond to the generalization of the BdG equations [25] to 
the case of spin-orbit coupled Fermi gases on optical lattices. 



C. Circulating mass currents 

Once we obtain the self-consistent solutions for the quasi- 
particle wave functions and the corresponding energy spec- 
trum, it is straightforward to calculate any of the desired ob- 
servables. For instance, in this paper we are interested in the 
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quantum mechanical probability-current (mass- or particle- 
current) of a fermions at site i defined by J a i = J^x + J^y, 
where Jj^ = —it l l ^ re {c\ i c r ,,i+e—H.c.) is the eth component. 
Using the Bogoliubov-Valatin transformation, we obtain 



J% 



-it. 



n,i-\-e 



n,z+e 



/(-e„) - H.c. 



(5) 



where we set the lattice spacing and h to unity. In our two- 
dimensional system, we expect mass currents to circulate 
around the central site, so that = J CT j0, where is the 
azimuthal angle. In addition, the local angular-momentum as- 
sociated with such a circulating mass current is in the z direc- 
tion, and its magnitude at a particular site i is simply related to 
the local mass current by L z ai = xi — y; J^, where Xi and 
yi are the coordinates of site i with respect to the center, i.e. 
minimum of the trapping potential. In this paper, we are inter- 
ested in the total absolute angular momentum per a particle 
which is given by L z a = (1/N a ) J2i n\Jai\- 

Having established the BdG formalism, now we are ready 
to present our numerical solutions for the ground-state and 
finite-temperature phases, which are obtained by solving 
Eqs. (f2}|4|i in a self-consistent manner. 



III. NUMERICAL RESULTS 

Our numerical calculations are performed on a two- 
dimensional square lattice with L = 48 lattice sites in each 
direction. We take Vq « 0.014i as the strength of the trap- 
ping potential, and discuss both symmetric (a = f3) and 
asymmetric (a ^ ft) SOC fields. Since we are mainly inter- 
ested in the low-filling population-balanced systems, we set 
Nf = Ni = 125, and study the resultant phases as a function 
of g, a and f3. However, we briefly comment on the effects 
of high filling and population imbalance on the ground-state 
phases toward the end of this section. 

In Fig. [U we show the strength of the j" mass current J^i at 
zero temperature (T = 0) for noninteracting (g = 0) systems 
with symmetric and asymmetric SOC, where we set a = ir/4 
and vary j3 parameter. We use Fig. [If a) as our reference data 
in all of our discussions below. First of all, the coarse-grained 
data near the trap center is a finite-size lattice effect and it 
does not affect our main findings. In addition, due to the 
time-reversal symmetry, J^i has exactly the opposite circu- 
lation for population-balanced systems and this current is not 
shown throughout this paper. As j3 decreases from the sym- 
metric case with /3 = 7r/4to/3 = 0, the asymmetry between x 
and y directions increase in such a way that the spin-flipping 
hopping gradually decreases to zero in the y direction, and the 
SOC field eventually becomes purely one dimensional in the 
x direction. We find that the mass current flows as long as 
the SOC field has both x and y components, and its strength 
decreases gradually as a function of increasing the asymmetry 
of SOC fields. Note that since the asymmetric SOC breaks the 
C4 symmetry of the square lattice, the resultant mass currents 
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FIG. 1: (Color online) The strength of the f mass current J-fi (in 
units of t) is shown as a function of lattice coordinates Xi and yi at 
zero temperature (T = 0) for noninteracting systems (g = 0). Here, 
we consider (a) a symmetric SOC where we set a = /3 — rr/4, 
and (b-c) asymmetric SOCs where we set ft to n/8 in (b) and n/16 
in (c), while keeping a = n/4. The lighter yellow (darker red) 
indicates clockwise (counterclockwise) circulation. J^i has exactly 
the opposite circulation and it is not shown. 



also have reduced symmetry in Figs. |TJb) and (c). In addi- 
tion, we see that the gas expands a little bit with increasing 
asymmetry, which is due to the decrease in density of single- 
particle states. 

In Fig. |2l we show the mass current and filling nfi of 
t fermions at T = for noninteracting systems with vary- 
ing symmetric SOC strengths a = (3. We find that first 
increases and then decreases as a — (3 is increased from 
to 7r/2. It is expected that the strength of the mass current 
rapidly increases as a function of increasing spin-flipping hop- 
ping with respect to the spin-conserving ones, since the pres- 
ence of a nonzero SOC is what allowed the mass current to 
form in the first place. However, as the SOC terms domi- 
nate, the chemical potential gradually drops below the band 
minimum, and then the system gradually crosses over to the 
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FIG. 2: (Color online) (a) The strength of the f mass current (in 
units of t), and (b) the filling of f fermions n-j-j are shown at T = 
as a function of x%, when yi — 0. Here, the system is noninteract- 
ing, and we set a = /3 to 7r/4 (solid black), 57r/16 (dashed red) 
and 77r/16 (dotted blue). Coarse-grained data is a finite-size lattice 
effect. 



molecular BEC side. We note that in contrast to the no-SOC 
case where the formation of a two-body bound state between 
t and I fermions requires a finite interaction threshold in a 
two-dimensional lattice, here it can form even for arbitrarily 
small interactions simply by increasing the SOC, due to the 
increased single-particle density of states. Therefore, it is ex- 
pected that the rapid increase in the strength of the mass cur- 
rent is followed by a gradual decrease at larger SOC, and that 
the mass current eventually vanishes at sufficiently large SOC. 
In addition, we see that the gas shrinks with increasing SOC, 
which is again due to the increased density of states. Note that 
the ratio between the trapping potential and the effective hop- 
ping terms also increases as a — (3 increases from to ir/2. 
We also calculate the absolute angular momentum per particle 
L z a and find approximately 0.53, 0.56 and 0.38 (in units of K), 
when a = f3 is set to 7r/4, 57t/16 and 7n/16, respectively. 
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FIG. 3: (Color online) (a) The strength of the f mass current J-j-j (in 
units of t), and (b) the filling of f fermions rifi are shown at T = 
as a function of Xi, when yi = 0. Here, the SOC is symmetric with 
a = ft = it/ 4, and we set g to (solid black) lAt (dashed red) and 
4.2t (dotted blue). Coarse-grained data is a finite-size lattice effect. 

In Fig. f3j we show the mass current JL and filling of t 



fermions at T = for symmetric (a = j3 = 7r/4) SOC with 
varying interaction strength g. The corresponding superfluid 
order parameter |A,| has an approximately peak value of 0, 
0.05£ and 1.4i at the trap center when g equals to 0, lAt and 
4.2i, respectively (not shown). It is clearly seen that increas- 
ing the g has qualitatively the same effect on the ground-state 
phases as increasing the SOC. For instance, we find that the 
mass current first increases and then decreases as a function 
of increasing g. In addition, the gas shrinks with increasing g, 
due to the formation of more tightly-bound Cooper pairs. We 
again calculate the absolute angular momentum per particle 
L% and find approximately 0.53, 0.49 and 0.22 (in units of h), 
when g is set to 0, lAt and 4.2i, respectively. Therefore, L z a 
again deviates significantly from w 0.5 as a result of increased 
filling around the trap center. 




FIG. 4: (Color online) (a) The strength of the f mass current J-j-j (in 
units of t), and (b) the filling of f fermions rifi are shown at finite T 
as a function of Xi, when yt = 0. Here, the system is noninteracting 
and SOC is symmetric with a — j3 — 7r/4, and we set T to (solid 
black) 0.14t (dashed red) and 0.35i (dotted blue). Coarse-grained 
data is a finite-size lattice effect. 

In Fig. [4] we show the mass current and filling of f 
fermions for noninteracting systems with varying temperature 
T, where we set symmetric (a = /3 = tt/4) SOC. Since finite 
temperature excites more and more particles to higher oscil- 
lator states as a function of increasing T, and this naturally 
leads to an increase in the system size, and decrease in the 
strength of the mass currents. This is clearly seen in the fig- 
ure, where the peak value of Ji^ reduces approximately to half 
of its T = value when T ~ 0.14i, but with a larger width. 
We calculate the absolute angular momentum per particle L z a 
and find approximately 0.53, 0.33 and 0.16 (in units of H), 
when T is set to 0, 0.14i and 0.35i, respectively. Therefore, 
L% decreases significantly from its ground-state value ss 0.5 
as the temperature increases. 

So far we only discussed population-balanced systems with 
a low filling at the trap center, i.e. n a i < 0.4, and here we 
briefly comment on the effects of high filling and popula- 
tion imbalance on the ground-state phases 112511 . First of all, 
due to the particle-hole symmetry of the Hamiltonian around 
n a i = 0.5, in addition to the mass-current peak circulating 
in the vicinity of the trap edge, an additional peak that is cir- 
culating mostly near the trap center is induced when the trap 
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center is close to a band insulator, i.e. n a i — 1. Therefore, 
mass currents show a double-ring structure in high-filling lat- 
tice systems. Second, the strengths of mass currents are quite 
robust against the effects of imbalanced populations. For in- 
stance, when the total number of J, fermions is reduced from 
N i = 125 in Fig. da) to N± = 20 while keeping N t = 125 
fixed, so that n-f « 0.3 and n± ~ 0.05 at the trap center, the 
maximum currents > J± « 5 x 10~ 3 t occur about the 
same radial distance as the population-balanced case shown 
in Fig. Ola). 



IV. DISCUSSION 

Having established the qualitative behavior of sponta- 
neously induced counterfiowing mass currents in spin-orbit 
coupled Fermi gases, next we argue that the origin of these 
currents can be understood via a direct correspondence with 
the chiral p-wave superfluids/superconductors. 



A. Correspondence with chiral p-wave systems 

The origin of spontaneously induced counterfiowing mass 
currents near the trap edge can be understood via a direct cor- 
respondence with the p x + ip a -superfluids/superconductors, 
e.g. A phase of liquid 3 He, for which it is known that a spon- 
taneous mass current also flows near the boundary |3,|5|]. In 
these p-wave systems, the mass current and hence the macro- 
scopic angular momentum is associated with the chirality of 
Cooper pairs. The chirality can be most easily seen by noting 
that the chiral p-wave order parameter At oc (x ± iy) ■ k, 
where k is the relative momentum of a Cooper pair, is an 
eigenfunction of the orbital angular momentum with eigen- 
value ±h. This mechanism explains our findings since it can 
be shown that the order parameter of spin-orbit coupled Fermi 
gases with Rashba-type SOC and s-wave contact interactions 
has chiral p-wave symmetry. To see this one needs to trans- 
form the Hamiltonian to the helicity basis, i.e. the pseudo- 
spin is no longer a good quantum number in the presence of 
a SOC, and the spin direction is either parallel or anti-parallel 
to the direction of in-plane momentum in the ± helicity basis. 
In the helicity-basis representation, it becomes clear that only 
intraband Cooper pairing occurs between fermions with same 
helicity, and that the order parameter for ± helicity pairing has 
Px T Wy symmetry 

Although the basic mechanism in chiral p-wave systems 
and spin-orbit coupled Fermi gases share some similarities 
with respect to spontaneously induced edge currents, they also 
differ in fundamental ways. For instance, in contrast to the 
chiral p-wave systems where the formation of Cooper pairs 
and hence the mass current requires an interacting system, 
in spin-orbit coupled Fermi gases currents flow even in the 
absence of interactions. In addition, the chiral p-wave sys- 
tems belong to the topological class of integer quantum Hall 
systems, and since these systems both break time-reversal 
symmetry, they exhibit spontaneous edge currents circulat- 
ing along a particular direction. However, spin-orbit coupled 



Fermi gases preserve time-reversal symmetry just like quan- 
tum spin Hall systems, and therefore, they both exhibit sponta- 
neously induced counterfiowing t and J. edge currents. Next, 
we comment on the experimental realization of spontaneously 
induced edge currents in atomic systems. 



B. Experimental realization 

In comparison to bulk (center) effects, there is no doubt that 
the edge effects are more difficult to observe in trapped atomic 
systems. This is because as the local density decreases from 
the central region to the edges, it gets harder to probe local 
properties. In addition, the local critical superfluid tempera- 
ture eventually drops below the temperature of the system near 
the edges, and the superfluidity is also lost starting from the 
edges inwards. One way to partially overcome such problems 
could be to load the optical lattice potential with higher parti- 
cle fillings such that the trap center is close to a band insulator, 
i.e. unit filling. In this case, an additional mass-current peak 
circulates mostly near the trap center (see above), and such a 
local current could be easier to probe with current experimen- 
tal capabilities. 

The measurement of the angular momentum of rotating 
atomic systems have so far only been achieved indirectly, by 
observing the shift of the radial quadrupole modes. While this 
technique was initially used forrotating atomic BECs 1261, 12711 . 
it has recently been applied to the rotating fermionic super- 
fluids in the BCS-BEC crossover [28]. We believe a simi- 
lar technique could be used for measuring the intrinsic angu- 
lar momentum of spin-orbit coupled Fermi gases, which may 
provide an indirect evidence for counterfiowing mass currents. 
In addition, we note that a realistic scheme has recently been 
proposed to detect topological edge states in an optical lattice 
under a synthetic magnetic field II29I1 . This proposal is based 
on a generalization of Bragg spectroscopy that is sensitive to 
angular momentum. 



V. CONCLUSIONS 

To conclude, we showed that the presence of a Rashba-type 
SOC spontaneously induces counterfiowing f and J. mass cur- 
rents. While these currents have a single peak that is circu- 
lating mostly near the trap edge in low-filling lattice systems, 
an additional peak that is circulating mostly near the trap cen- 
ter is also induced in high-filling lattice systems, exhibiting a 
double-ring structure. We note that our results for the low- 
filling lattice systems are directly applicable to the dilute sys- 
tems (without the optical lattice potential), for which we ex- 
pect qualitatively similar behavior. These currents flow even 
in noninteracting systems, and they are quite robust against 
the effects of imbalanced populations and/or asymmetric SOC 
in x and y directions. However, their strength decreases to- 
ward the molecular BEC limit, which can be achieved either 
by increasing the SOC or the interaction strength. We argued 
that the origin of spontaneously induced counterfiowing mass 
currents in spin-orbit coupled Fermi gases can be understood 
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via a direct correspondence with the chiral p-wave superflu- 
ids/superconductors, however, with some fundamental differ- 
ences in between. 
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